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Abstract 


In the analysis of the vibration response and structure-borne 
vibration transmission between elements of a complex structure, statistical 
energy analysis ( SEA) or finite element analysis (FEA) are generally used. 
However, an alternative method is using vibrational power flow techniques 
which can be especially useful in the mid frequencies between the optimum 
frequency regimes for FEA and SEA. Power flow analysis has in general been 
used on one-dimensional beara-like structures or between structures with 
point joints. In this paper, the power flow technique is extended to 
two-dimensional plate like structures joined along a common edge without 
frequency or spatial averaging the results, such that the resonant response 
of the structure is determined. The power flow results are compared to 
results obtained using FEA at low frequencies and SEA at high frequencies. 
The agreement with FEA results is good but the power flow technique has an 
improved computational efficiency. Compared to the SEA results the power 
flow results show a closer representation of the actual response of the 
structure . 


Introduction 


The use of power flow techniques for the analysis of the vibration 
response and structure-borne vibration transmission between elements of a 
structure has in general been restricted to either one-dimensional beam 
type structures [1,2] or to structural elements with only point attachments 
[3]. The problem of line joints between structural elements has only been 
studied using power flow techniques averaged over defined frequency bands 
(one-third octave frequency bands) [4] in which case the results that are 
obtained are identical to those of statistical energy analysis (SEA). 

Thus, the potential of the structural power flow method to give results for 
the detailed structural response for the elements of the structure, which 
as shown in [l] can be very important in the mid frequency range is 
completely lost. 

The basic concepts of the power flow method and its application to 
determine the transmitted power have been demonstrated in a number of 
published papers [1,2]. In the previous paper by this author comparisons 
were also made between the results obtained using the power flow technique 
and results obtained using both SEA and finite element analysis (FEA). As 
shown in that paper [1] while FEA and SEA are excellent tools to be used in 
the low and high frequency regimes, respectively, a gap is left in the 
mid-frequency range where FEA becomes computationally unwieldy, and SEA 
would only give mean response levels with the actual response of the 
structure, especially near resonances, being potentially significantly 
different from these mean levels. The gap in the mid-frequency range, 
between FEA and SEA, can be satisfied using the power flow techniques. In 
fact since FEA and SEA are well established tools, these can be used to 
determine the results at low and high frequencies. These results are then 
compared to results obtained using the power flow method in these frequency 
regimes and by ensuring that there is a good match between the two sets of 
results in each respective frequency range, this can serve as a 
verification of the power flow model. 

Following is a brief general outline of the power flow method. The 
structure is modeled by a series of coupled substructures with forces or 
moments introduced at the joints between the substructures. The power flow 
into the source substructure and between the substructures is expressed in 
terms of input (point, line or area) and transfer structural mobility 
functions. The definition of structural mobility is the ratio of rate of 
charge of displacement (translational or rotational) per unit input load. 
For point mobilities the response and the load are measured at the same 
location, while for transfer mobilities the response and the load are at 
separate locations. The exact definition of the mobility functions to be 
used at the different joints between the substructures will depend on the 
type and configuration of the joint. The detailed response for the 
structural response of the elements of the structure is obtained by 
ensuring that the mobility function used in the analysis are not spatially 
or frequency averaged. 


The use of the power flow analysis for the determination of the 
response of a two-dimensional plate-like structure is presented in this 
paper. In the analysis two plates are joined at right angles along a common 
edge and one of the plates, the source plate substructure is excited by a 
point force. The transmitted vibrational power, to the attached second 
plate, the receiver plate, presented both per unit input force and as a 
ratio to the input power, is computed. The results obtained using the 
power flow technique are then compared to results obtained using FEA and 
SEA. However, in this case, unlike what was done in [1] the comparisons are 
restricted to low and high frequency regimes respectively. With FEA the 
computation time required to obtain the transmitted power as a function of 
frequency were becoming prohibitively long as the frequency range 
increased. 


L-Shaped Plate Structure 

Using the power flow method, the vibrational power transmitted between 
two plates joined along a common edge to form an L-shaped plate is 
computed. One of the plates is excited by a point load arbitrarily 
located, but initially the excitation location is kept away from the common 
edge. The L-shaped plate, coordinate system and loading are shown in 
figure (1). The following conditions were assumed for the two plates. 

1. The plates are thin compared to the wavelength so that 
inplane forces, rotary inertia and shear can be neglected. 

2. The joint between the two plates is a pinned joint. 

These conditions were only assumed such that some boundary conditions can 
be defined- However, other boundary conditions and the analysis for thick 
plates, where inplane forces and rotary inertia and shear are important can 
be analyzed in a similar fashion. In fact, thick plates with different 
boundary conditions are presently being investigated. 

Using the power flow method, a model for the structure is first 
defined. In general, a power flow model consist of a number of source and 
receiver substructures for which ways of obtaining the mobility functions 
are available. For the L-shaped plate structure considered here, the 
structure is divided into two substructures which are joined at the edge. 
For the L-shaped plate structure considered here, the structure is divided 
into two substructures which are joined at the edge. In general terms, 
using the two substructure model with one force input, shown in figure (2) 
the expressions for the input and transmitted power are [2] 


P. 

xn 


— Real 


M 12 M 21 

M2 + M3 


I F (f) | 2 


2 


( 1 ) 


2 

Real | M 3 j- 

(2) 

where M^, i - 1 , 2,3 are the input mobilities at locations 1 , 2 and 3 and 
M^2 and M21 are the transfer mobilities. 

In the case of the L-plate structure with point load excitation, 
would be the input mobility for the source plate substructure when isolated 
from the receiver substructure. The input mobility here would be the 
normal surface velocity at the excitation location per unit applied load. 

M2 and M3 represent input mobilities for thesource and receiver plate 
substructures respectively, at the location of the joint between the 
substructures. However, because the plates are joined along an edge, the 
definition of these mobility functions is different from that of . M2 and 
M3 will represent the rate of change of displacement, translational and 
rotational at all the points along the junction edge of the plate, per unit 
applied inplane, transverse and moment loads. The motion at every location 
along the junction will be dependent not only on the loading at that 
location but on the load distribution along the entire edge. To ensure that 
the results to be obtained from the analysis are relevant to the solution 
to the response of the global structure, the distribution of .the loading 
along the edge must be identified. M^2 and M21 represent the transfer 
mobilities between the excitation location and the joint location 
respectively and vice versa. In this case as well, similar to M2 and M3, 
because the joint is along one edge the transfer mobilities will be defined 
as follows. M12 is given by the ratio of the rate of change of 
displacement (also translational and rotational.) at every location along 
the joint and the applied loading. Thus M32 will in general be a function 
of the location along the edge. M21 is given by the ratio of the rate of 
change of displacement at the location where the excitation is applied per 
unit applied loading along the edge. 

Since some of the mobility functions described in the previous 
paragraphs are function of the location along the edge, the power 
transmitted will thus be given by an integral over the entire length of the 
joint . 

Thus far, the discussion on the required mobilities, has been 
restricted to general definitions in the context of two-dimensional 
structures joined along a common edge. Since the assumptions are being 
made that the plates are thin compared to the wavelength and that no 
transverse motion is allowed at the joint, some simplifications can be made 
to the definition of the mobility functions. For thin plates the main type 
of motion under transverse loading is bending or the generation of flexural 
waves. If only transverse displacements are considered, the excitation 
along the joint edge will be in the form of a distributed moment. Thus the 
parameters that one is dealing with for the particular case studied here 
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are, a point load, transverse displacement at the point of application of 
this load, a distributed moment along one edge and rotational displacement 
along this same plate edge. 

Having developed the power flow model for the problem and identified 
the required mobility functions and their respective definitions, the next 
step is to develop expressions for these mobility function. In developing 
these mobility expressions, it is not necessary to consider the global 
structure but each substructure is considered separately. For the L-plate 
problem, each substructure consists of a flat plate structure and thus in 
developing the mobility functions a flat plate with the different types of 
loading as described above needs to be considered. Thus the analysis of 
the L-plate structure is simplified to the analysis of flat plate 
structures . 

In deriving an expression for , a plate with the loading as shown in 
figure (3(a)) is considered. The boundary conditions of the plate are the 
same as those for the global structure, that is all edges are pinned or 
simply supported and therefore the boundary condition are zero displacement 
and zero bending moment. Since a thin plate is being considered, the 
response can be expressed in terms of a series solution with assumed mode 
functions that satisfy the defined boundary conditions. Thus for [5} 
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where f is the frequency, a and b are the width and length of the plate, p 
is the plate material density and h is the plate thickness. f n is the nth 
damped natural frequency of the plate given by 


where D* is the plate complex flexural rigidity which includes structural 
damping and n represents the double subscript (n^, n2 ) • 

An expression for M]^ can also be derived from this plate 
configuration since M^2 represents the angular displacement along one edge 
of the plate, which would be the joint edge per unit applied force when a 
transverse load is applied somewhere on the plate. That is 
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To derive expressions for M2, M3 and M21 a plate with the 
configuration shown in figure ( 3 (b)) is considered. In this case a moment 
is applied along on edge of the plate. The boundary condition are 
therefore simple supports along three edges, zero displacement and zero 
moment, and one edge pinned with an applied moment. Because of these 
uncommon boundary conditions the first task is to derive an expression for 
the mode functions and then use a series solution to derive the mobility 
functions. The displacement u (x,y) can be given by an expression of the 
form [ 6 , 7 ] 
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where Uj ( y ) is a function of the boundary conditions in the y-axis 
direction. This function must satisfy the conditions 
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where T (x) is the applied moment along the plate edge. Using these 
boundary conditions and the general solution for the response in one 
direction, one can derive an expression for (y> to give, (the derivation 
of the following equations can be found in Appendix A) 
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and the wavenumbers k^ and k^ are defined by [6,7] 
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where k x and ky are the plate flexural wavenumbers in the x and y 
directions respectively. Thus k x and k y satisfy the relationships 
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All wavenumber components are complex quantities to include the structural 
damping. 


Equation (8) can be used to obtain M2, M3 and M2 1 - However, before 
any of these mobility functions can be evaluated, an expression for T n is 
required, which would require knowledge of the distribution of the moment 
along the plate's edge. 

The distribution of this moment is a function of the type and location 
of the loading of the global structure. One way around this moment 
distribution dependency is to derive an overall function for 
M12/ (M2 + M3) to be used in the power flow expression, equation (1) and 
then a separate expression for #21 m Also, each of these functions would 
be derived for each mode separately and the overall result obtained by 
summing all the contributions from all the modes. Thus after some 
manipulations starting from equation (8) the following relationships can be 
derived. 
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where the kj and k2 in the equation for A for the source substructure, are 
different from those in the equation for B which is for the receiver 
substructure, and 







sin(k„y ) 
2 o 


sinh (k y ) 


, j / „ n, 7TX 

^7] s Kt^) 


aVphD* *-sin(k 2 b 1 ) sinh (k 


( 13 ) 



and for the particular case of two identical plates the overall power input 
obtained by substituting equations (3), (12) and (13) into equation (1) is 
given by 
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To obtain an expression for the transmitted vibrational power it is 
more important that the distribution of the moment along the joint edge is 
known. The transmitted power is given by the integral of the transmitted 
power per unit length over the entire length of the joint. For the special 
case of two identical plate substructures of equal size the expression for 
the transmitted power is given by 
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Both expressions for the input and the transferred power are functions of 
the input loading position (x Q , y Q ) and the material charac teristics of the 
plates given by the density, thickness, area and natural frequencies. 


Results 


The input and transmitted vibrational power for an L-shape plate 
structure were computed using the above derived equations. However, it is 
necessary to verify the results and thus comparison is made with results 
obtained using FEA and SEA. Before proceeding to discuss the results a 
brief description of the FEA model is given. The package used for the FEA 
is the MARC finite element package [8]. 

The model for the L-plate consists of 63 nodes and 48 flat thin plate 
elements (element type 50 in MARC). The decision on the number of nodes 
was made after computing the required number of nodes required to achieve 
good accuracy in estimating the mean square velocity over the entire 
surface of the structure, up to a maximum frequency of 1 KHz. The total 
number of nodes was however slightly decreased to reduce the computation 
time. With 63 nodes good accuracy is expected up to 600 HZ. In fact this 
is the role to be played by FEA as a verification of the results at low 
frequencies. The FEA model of the plate is shown in figure (4). This is 
the model used for the modal analysis which had two extra rows of elements 
to obtain a better indication of the behavior (rotation) of the joint. 
During modal analysis the computation time is rather short and thus the 
addition of the extra two rows of elements were not a problem with computer 
run time. 

FEA were used to compute both the natural frequencies and mode shapes 
and the response of the receiver plate when a harmonic force is applied to 
the source plate. Since this analysis was intended to be a verification of 
the power flow results, the same plate structural characteristics and point 
of excitation were used. The modal analysis was performed to ensure good 
correlation between the FEA model and power flow model. 

In the harmonic analysis the model with 63 nodes was used. The 
objective of this multi-increment frequency analysis was to compute surface 
deflection data from which the spatial average velocity squared is 
obtained. The power transferred to the receiver plate, which is equal to 
the dissipated power by conservation of power, is given by 

p trans " n2irf pabh< | V2(f ) | > (16) 



where rj is the receiver plate loss factor and < > denotes spatial 
averaging. In the spatial averaging of the velocity squared of the 
receiver plate, only the internal nodes were considered. 

A number of excitation location were considered and two typical 
results for the transmitted power per unit force squared are shown in 
figures (5(a) and (b)). On the same figures are shown the power flow 
results and the results using the FEA. The frequency resolution used for 
the FEA is different for figures (5(a)) and (5(b)). For figure (5(a)) a 
resolution of 1 Hz was used in the frequency range between 1Hz and 100Hz 
and 10Hz resolution for the frequency range 110 to 1000Hz. In figure 
(5(b)) the frequency resolution is constant at 5Hz . 

The results obtained using the power flow method for the transferred 
power were also compared to results obtained using SEA. With SEA a ratio 
of the transferred power to input power is usually given thus a ratio 
between the transferred and input power is obtained by dividing equation 
(15) by equation (14). Two sets of results for excitation in the middle of 
the source plate are shown in figure (6). The power ratio curve computed 
using the power flow technique show some anomalies, in that the ratio 
exceeds unity. These are due to resolution errors and infact occur close 
to the resonant peaks, where the exact definition of the peak is important. 
Also the computational method is breaking down at two frequencies, just 
above 400Hz and 800Hz . Both of these errors are currently being 
investigated. 

As can be observed in figures (5) and (6) , good agreement is 
obtained between both the power flow and the SEA and FEA results. However, 
the computation for power flow is much more efficient compared to FEA when 
the spectrum of the transmitted power is required. Compared to SEA power 
flow is computationally more intensive, but the results obtained are closer 
to the actual response of the structure especially at the resonant 
frequencies . 


Conclusion 


In conclusion, a power flow technique has been presented in this paper 
which allows the computation of the response of a complex structure. What 
has been shown in this paper is the extension of the power flow technique 
to two-dimensional structures, and the results obtained match very closely 
with results obtained using FEA, but with potential saving in computation 
time. Also, the results are compared to results obtained using SEA and 
while the latter only predicts mean values for the response, the power flow 
method gives a better representation of the actual response of the 
structure . 

In this paper, only the formulation of the power flow method for 
two-dimensional structures have been presented. However, having shown that 
the power flow method is a very powerful technique, and that it is not 
limited to one dimensional or solely point connected structures but can be 


used for the analysis of two-dimensional structure, the influence of 
parameters of the structures and the excitation on the response of 
individual substructures will be investigated. Because of the fact that 
the power flow results represent the narrow band spectrum for the response, 
the effect of the structure and excitation parameters on the exact response 
of the structure can be deduced. 

Additionally, the technique needs to be extended to more complex 
structures. This would require that the respective mobility functions used 
in the analysis be better defined than they are here. The shape of the 
moment function at the joint needs to be fully defined from which mobility 
functions for various condition can be derived and stored in a data base 
for future use. The analysis of a complex structure would then be an 
exercise in putting together the right mobility functions which would be 
available in the data base. The choice of the mobilities can be verified 
by comparing the very low frequency results and modal frequencies to 
results obtained using FEA, thus verifying the power flow model and the 
mobility functions used. 
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L-PLATE 



Thickness: 0.00635 m 

Dimensions: 1.0 m x 0.5 m 

Loss Factor: 0.01 

Figure 1. Plate structure showing plate characteristics. 



Figure 2. General power flow model. 


simply supported 



Figure 3. Power flow model substructural components - of the 
L-shaped plate. (a) Plate with point excitation; 
(b) Plate with distributed edge moment excitation. 
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Figure 5. Results of the power flow analysis and the FEA for the 
transferred power, (a) Excitation in the center of the 
source plate; (b) Excitation off center. 
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plate, (a) Comparison of transferred and input power; 
(b) power ratio and SEA results. 











Appendix A 


The mobility functions for the case of a plate simply supported along 
three edges, and pinned with a distributed moment applied along the fourth 
edge are derived as follows. 

The displacement normal to the surface of the plate is given by 
equation (6) which satisfies the boundary conditions (see figure 3(b)), 


u(a,y) = 0 


The function *Ky) in equation (6) is chosen such that the conditions 
given by equation (7) are satisfied. The expression for \|)(y) is obtained 
by solving the one-dimensional differential equation (V4 - k4)u(x,y)=0 
in the y-axis direction. That is [7] 

IKy) * A cosMkqy) + B sinMkjy) + C cosO^y) + D sin(k 2 y) (A. 2) 

where A, B, C, D are constants dependent on the boundary conditions in the 
y-axis directions, and and k£ are defined in equations (10) and 
( 11 ). 
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From the boundary conditions, 

A = C 5 0 (A. 3) 

and to solve for B and D the distributed moment is spatial transformed into 
its wavenumber (or wavelength) components 
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from which equation (8) is derived. 

The distributed moment dependency from the mobility functions is 
eliminated by taking ratios of the mobilities. From the basic equation 
relating the motion of a coupled structure 

02<x,b) = e F (x,b) + M2 (x) T(x) (A. 7 ) 

where 0 2 (x,b) is the angular velocity at the junction; 

0 F (x,b) is the angular velocity at the junction due to the 
application of the force F; 

and M 2 (x) T(x) is the angular velocity due to the applied moment. 

But the angular velocity at the junction roust be the same for both 
substructural components, 

02(x) = 03(x) « -M 3 (x) T(x) (A. 8) 


and thus 


0 F - -(M 2 + M 3 ) T(x) (A. 9 ) 

where the x dependency of M 2 and M3 is suppressed for brevity. For 
M 2 and M 3 
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and similarly for M3, where for M2 i*l and for M3 i* 2 . Along the 
x-axis direction both substructures have the same mode number. Thus 
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and from equation (A. 9) 
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The advantage of £(y) over ty(y) is that £(y) is independent of 
the distributed moment. From equation (A. 12) 
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Substituting in equation (A. 14) for the numerator by equation (5) and for 
the denominator by equation (A. 13), equation (12) is obtained. 


From definition 


M 21 (x) * V f /[T(x)] 


(A. 15) 


where Vp is the velocity at the excitation point due to the moment 
applied along the edge. Thus from equation (8) 
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By combining equations (A. 14) and (A. 16) the expression for the input 
power (equation 14) is obtained. 



The transmitted power to the receiver substructure is given by the 
integral over the entire length of the joint: 
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For the case of two identical plate substructures, using the basic 
relationship 



(A. 19) 


Both terms within the integral are sinusoidal functions and represent 
different modes of the structure for different values of nj . Due to 
orthogonality of the modes equation (A. 19) simplifies to: 
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Substituting in equation (A, 20) by equations (A. 13) and (8), the 
expression for the transmitted power, equation (15) is obtained. 



Appendix B 


Exp e r ime n tal Analysis Of Power Transmission In An L-shaped Beam 


Introduction 


An experimental analysis was performed on an L-shaped beam structure 
(figure B.l). The objective of the experiment was to measure the power 
transmitted from the source beam to the receiver beam, which was excited 
by a harmonic force at one end. The results are compared to those obtained 
in the analytical study [1]. The experimental procedure is to measure the 
response at several points on the receiver beam and to obtain a spatial 
average of the square of the normal surface velocity over these points. In 
using this method, it is assumed that the power transmitted to the 
receiver beam is equal to the power dissipated since the beam is not 
connected to anything other than the source beam. The time averaged power 
dissipation in a finite structure is given by equation (16). 

Experimental Setup 

The experimental setup consists of two identical steel beams of 
thickness 3/8 inch (0.0095m), width 3 inches(0. 0762m) , and length 36 
inches (0.9144m). The beams are welded to a 5/8 inch (0.0159m) diameter 
steel rod such that these form a 90 degree angle at the joint (figure B.l) 
The rod provides a means of pinning the joint so that the experimental 
setup matches the boundary conditions used in the theoretical analysis. 

The rod was held by two flange bearings mounted to a base structure. The 
harmonic excitation is provided by means of a Wilcoxon electromagnetic 
shaker attached to the source beam via an impedance head and supported by a 
heavily damped support rig. 

The experiment consists of two parts: 

1) Determination of the damping loss factor 

2) Measurement of the spatial average velocity of the beams and 
the processing to obtain the dissipated power. 

The structural damping was measured using both the decay method and 
the half-power bandwidth method. The measurements were performed using a 
calibrated impact hammer, with an accelerometer at the end of the second 
beam. This measurement point was chosen because it does not represent a 
node at any of the natural frequencies. From damping measurements on the 
undamped L-beam it was determined that the bearings and the exciter 
attachment to the beams were the main source of energy dissipation. Since 
this would create undesirable effects in doing the measurements, 
viscoelastic material was added to one side of each beam. 

Damping measurements were again performed for the damped L-beam when 
freely suspended and when mounted in its support rig and the shaker 
attached. For the latter case the damping was measured using a single 
frequency excitation and measuring the decay after the excitation is 


switched off. The results of the damping measurements are shown in Table 
(B.I). It was observed that the addition of the viscoelastic material 
increased the damping sufficiently such that the influence of the bearings 
had little effect on the structural damping. 

In the second part of the experiment, the transverse acceleration 
response due to band limited white noise excitation was measured at ten 
uniformly spaced points on the receiver beam. The acceleration measurements 
were made using a B&K 4375 accelerometer. The data was processed 
according to the following procedure. 

1) A transfer function was first obtained between the spectra of 
the acceleration response data and the force data using an HP 5451C 
structural dynamics analyzer. 

2) The frequency domain data were then divided by j 2 7rf to obtain the 
velocity response and and then squared* 

3) These two steps were performed at each of the ten measurement 
points with the resulting ten sets of data averaged to obtain the spatial 
averaged velocity squared per unit input force squared. The frequency 
range used was between 0 and 1.25 KHz. 

4) Finally the data was multiplied by the measured damping loss 
factor and the beam mass and frequency to obtain the dissipated power. 

Results 


The measurement results for the power transmitted to the receiver beam 
are compared with the results given by the power flow method in figure 
(B.2) There are several comments which can be made about these results. 
First there is good agreement between the measured data and the analytical 
predictions in terms of the natural frequencies of the coupled system. In 
general any offset of peaks in the measurement curve from the peaks in the 
power flow curve may be attributed to insufficient frequency resolution in 
the measurements. Some discrepancies in level are obtained at the 
resonant frequencies which also may be due to resolution. Below 100 Hz it 
is difficult to draw any conclusions about the measured data. The cutoff 
frequency of the electrodynamic shaker and power amplifier is 20 Hz which 
makes the data below this frequency unreliable. The discrepancies in the 
results in the 40 to 100 Hz range are attributed to vibration of the 
supporting structure, both for the beam and the electrodynamic shaker. 
Additionally there are two peaks (at 112 Hz and 218 Hz) in the measured 
data which do not appear in the analytical; results. These are attributed 
to support structure resonances. However in general there is reasonable 
agreement between the analytical and the experimental results. 

Work is still in progress to measure the power flow vectors along the 
two beam substructures. Two techniques are being used which differ in the 
number of accelerometers used in the measurements. Some results have been 
obtained with a two accelerometer arrangement [B.l], however the results 



obtained thus far have not been fully explained. A possible reason for the 
unexpected results is the very large standing wave ratio which is the ratio 
of the amplitude of the standing wave to the propagating wave. This may 
possibly not allow us to do these type of power flow measurements using 
this arrangement. An alternative is to use a three accelerometer 
configuration [B.2] however no results have yet been obtained using this 
technique . 
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TABLE B.l 


Measured Damping Loss Factors 


Mode // 

Free 

Undamped 

Free 

Damped 

Pinned 

Damped 

Pinned, damped 
(shaker Attached) 

i 

0.005 

0.011 

0.007 


2 

0.004 

0.012 

0.008 

0.008 

3 

0.003 

0.008 

0.006 

0.010 

4 


0.008 

0.009 

0.018 

5 


0.008 

0.008 

0.020 

8 


0.009 

0.012 

0.017 

9 


0.009 

0.012 

0.018 
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Figure B.l. L—shaped beam structure. 





